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We analyze the topological structure of quenched QCD in the presence of static color sources. Distributions 
of the topological charge density around static quarks and mesons are computed in both phases of QCD. We 
observe a suppression of topological fluctuations in the vicinity of external sources. In the confinement phase, the 
suppression occurs in the whole flux tube between the static quark and antiquark. 
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1. The physical mechanism of quark confine- 
ment in QCD is still not completely understood. 
It is generally believed, however, that it is due to 
non-trivial topological properties of the theory. 

To get some insight into the confinement mech- 
anism and the corresponding vacuum structure it 
is instructive to study various physical quantities 
locally around color sources. In previous works 
we investigated the chiral condensate, charge den- 
sity, color field distribution, etc., around static 
sources in full QCD at finite temperature [1]. We 
observed that the vacuum structure changes sig- 
nificantly around external charges. Our results 
agree with phenomenological models like the bag 
model where the vacuum inside the bag is of per- 
turbative nature and differs from the nonpertur- 
bative medium outside. 

In the present work we focus on the topologi- 
cal properties of QCD. We study topological fluc- 
tuations around static quarks and static mesons 
in both phases of quenched QCD by computing 
correlation functions between Polyakov loops and 
the topological charge density. 

2. There is no unique way of defining a lat- 
tice operator that converges in the naive contin- 
uum limit to the continuum expression for topo- 
logical charge. We have chosen two possibilities, 
the "plaquette" (P) and the "hypercube" (H) 
definition [2], in which the topological charges 
Q( P ' H ) are defined as integrals (sums) of topo- 



logical charge densities: 
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where operators O^'^J are (in usual notation) 

OiV P a = U^(x)U pa (x), (2) 

Ofilo = U tl (x)U v {x + p)U p (x + p + v) 

x U a (x + p + v + p)Ul(x + v + p + a) 

x ul(x + p + a)Ul(x + a)Ul(x). (3) 

At finite coupling [3 and finite lattice spac- 
ing neither of the above definitions yields inte- 
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Figure 1. Squared topological charge density 
around a static quark after 0, 5, and 10 cooling 
steps at [3 = 5.8 (deconfinement phase). 
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Topological Charge Density in a Static Meson 
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Figure 2. Squared topological charge density around a static quark-antiquark pair at [3 = 5.6 and [3 = 5.8: 
(a) d = 2, after and 5 cooling steps, (b) d = 4, after 5 cooling steps. 
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ger topological charges. This is known to be due 
to renormalization factors relating the continuum 
and lattice operators of topological charge [3]. 
This problem can be overcome by using a cooling 
procedure. Cooling is subsequent minimizing of 
the action with respect to all individual links. We 
employ a slightly modified cooling method from 
Ref. [4]. Within a cooling step, we minimize in 
all three diagonal SU(2) subgroups of SU(3), in 
random order, and finally check after any mini- 
mization step that the changes of links are small 
enough not to destroy immediately the topologi- 
cal charge. After a few cooling steps, topological 
charges of lattice configurations, determined us- 
ing both definitions, agree and are close to integer 
values. 

A static quark is represented by the Polyakov 
loop L(r) which describes the propagation of a 
charge with infinite mass. To obtain the local dis- 
tribution of the topological charge density around 
single quarks and around a meson with quark- 
antiquark separation d we calculate the correla- 
tion functions (L(0)q 2 (r)) and (L(O)L^ (d)q 2 (r)), 
respectively, and normalize them to the corre- 
sponding vacuum expectation values. We take 
the square of the topological charge density as the 
simplest choice leading to nontrivial correlations. 

3. The simulations were performed on an 8 3 x 4 
lattice with periodic boundary conditions using a 
Metropolis algorithm. We evaluated path inte- 
grals in pure SU (3) with standard Wilson action 
both in the confinement (/? = 5.6) and in the de- 
confinement phase (j3 = 5.8). We made 10 5 iter- 
ations per [3 and measured our observables after 
every 50th iteration. Each of these 2000 configu- 
rations was subjected to 50 cooling steps. 

In Fig. 1 the square of the topological charge 
density around a single static quark in the decon- 
finement phase is shown after 0,5, and 10 cooling 
steps, for both definitions of the density. Without 
cooling practically no effect is visible, but cooling 
soon uncovers a suppression of the density. Re- 
sults from both definitions of the density agree af- 
ter a few cooling steps. In the confinement phase 
one expects (L(0)q 2 (r)) = 0, what was confirmed 
by our measurements. This reflects the fact that 
a single quark is a forbidden state in the confine- 
ment phase of pure gauge QCD. 



In Fig. 2 distributions of the topological charge 
density (squared) around a meson are presented 
for both phases of pure SU (3). The data for and 
5 cooling steps are visualized at (/(/-distance d = 2 
in Fig. 2a. For 5 cooling steps the suppression be- 
comes again clearly visible. The results indicate 
a stronger effect in the confinement regime. The 
difference between both phases is more striking 
in Fig. 2b where plots for the 5th cooling step 
and larger separation, d = 4, are drawn. The val- 
ley of suppressed topological charge density in the 
confinement phase reflects the formation of a flux 
tube between the static quark-antiquark pair. In 
the deconfinement phase, no flux tube is visible 
and the effect is peaked near the position of the 
sources (cf. Fig. 1). 

4. Our results show a change of the topological 
properties of lattice vacuum configurations across 
the deconfining phase transition. Configurations 
with nonvanishing topological charge appear very 
seldom in the deconfinement phase. The topolog- 
ical charge density is lowered in the vicinity of 
external color sources in both phases of quenched 
QCD. However, the effect is more pronounced 
in the confinement phase where it is observed in 
the whole flux tube between the static quark and 
antiquark. A similar investigation for QCD with 
dynamical fermions is in progress. 
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